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Abstract 



In the dilute instanton gas model of the QCD vacuum, one expects a strong 
spatial correlation between chirality and the maxima of the Dirac eigenvectors 
with small eigenvalues. Following Horvath, et al. we examine this question 
using lattice gauge theory within the quenched approximation. We extend 
the work of those authors by using weaker coupling, j3 = 6.0, larger lattices, 
IG'^, and an improved fermion formulation, domain wall fermions. In contrast 
with this earlier work, we find a striking correlation between the magnitude 
of the chirality density, {x)'^^iIj{x)\, and the normal density, ijj'^ {x)iIj{x), for 
the low-lying Dirac eigenvectors. 
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I. INTRODUCTION 



In a recent paper, Horvath, Isgur, McCune and Thacker suggest that an important test 
of various theoretical models of the QCD vacuum can be made by examining the degree to 
which the space-time localization of a low-lying eigenmode of the Dirac operator is correlated 
with non-zero chirality. Such correlations are expected in the dilute instanton gas model of 
the QCD vacuum. 

In this model, gauge configurations composed of widely separated instantons and anti- 
instantons support low-lying Dirac eigenvectors which are near superpositions of the zero 
modes that these (anti-)instantons would posses if in isolation. These near zero modes 
are proposed |0| to provide the non-zero density of Dirac eigenvalues, p(A) for vanishing 
eigenvalue A, required by the Banks- Casher formula and the non-zero QCD chiral con- 
densate, 7rp(0) = -12(gg) = (245MeV)^ where we give the value for [qq) obtained in a 
recent quenched calculation |Q and the factor of 12 comes from our particular normalization 
conventions. Since the isolated zero modes associated with an (anti-)instanton are entirely 
(left-)right-handed, one should expect a strong correlation of handedness with the locations 
at which ipl^ipn{x) is large. 

Thus, this class of instanton-based models of the QCD vacuum can be tested by search- 
ing for such strong correlations between the chiral density ipl^'y^ipnix) and normal den- 
sity ipl^ipnix) for the low- lying eigenmodes of the Dirac operator on a configuration- by- 
configuration basis within a lattice QCD calculation. As emphasized by Horvath, et al. 
earlier arguments of Witten suggest that the large Nf. behavior of the rj' mass may be 
inconsistent with the predictions of such an instanton model, providing additional motivation 
for such qualitative tests of the instanton picture. 

In the above cited work of Horvath, et al., the authors search for such correlations using 
30 gauge configurations obtained in the quenched approximation on a 12^ x 24 lattice volume 
with Wilson fermions and the Wilson gauge action with /3 = 5.7. In this paper, we extend 
their study to smaller lattice spacing, = 1.922 GeV, using a 16^* lattice volume which 
has a similar physical size of ~ 1.6 Fermi. Making use of configurations that we had already 
analyzed, we discuss here results from two different gauge actions. The first is the standard 
Wilson action with /3 = 6.0 for which we have 32 configurations and the second uses the 
Iwasaki action with jS = 2.6, a value tuned to achieve the same lattice spacing. The latter 
ensemble contains 55 configurations. 

Finally we employ an improved fermionic action, using the domain wall formalism of 
Shamir |^. Further references to this method, our implementation of this lattice fermion 
action and a complete description of the notation used in this paper can be found in Ref . . 

As is detailed below, using this finer lattice spacing and improved fermion formulation, 
we find a very different result from that seen in the earlier work of Horvath, et al. The 
chirality of the n*'' eigenmode, \ipl^{x)'y^ilJn{x)\ is surprisingly close to maximum at those 
lattice sites x where the eigenmode is large. In Section |I| we briefly review the properties 
of the spectrum of the continuum Dirac operator and the closely related hermitian Dirac 
operator that we study. Section [TTTI explains the diagonalization procedure. In Section [TVI we 
present our results and analysis of the domain wall fermion eigenvectors. Finally, Section 
contains a brief conclusion. 
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II. PROPERTIES OF THE CONTINUUM DIRAC OPERATOR 



In the continuum, the Euchdean Dirac operator is usually written as a non-hermitian sum 
of an anti-hermitian operator Ip = 7^^/}^ where is the usual gauge covariant derivative 
and a real mass term. The combined operator, Ip + m, is easy to analyze being the sum 
of commuting hermitian and anti-hermitian pieces. However, the standard Wilson and 
domain wall lattice Dirac operators are more complex because the so-called Wilson mass 
term also contains derivatives causing the hermitian and anti-hermitian parts of these lattice 
Dirac operators to fail to commute. For the Wilson operator, Dw, this is easily remedied by 
working with the product which is readily seen to be hermitian. A similar construction 

is possible for domain wall fermions where the product 'j^Rt^D is hermitian. Here D is 
the domain wall Dirac operator, given for example in Eq. 1 of Ref. 0, and the reflection 
in the midplane in the fifth dimension, s — > — 1 — s. 

Our treatment differs from that of Horvath, et al. in that we work with the hermitian 
operator, finding real eigenvalues and their corresponding eigenvectors while the authors of 
Ref. [|I| determine the complex eigenvalues of the non-hermitian Wilson operator and their 
corresponding eigenvectors. In order to compare our two results, we now review the relation 
between the eigenvectors of corresponding continuum operators: Ip -\- m and 7^(-^ + m). 

We begin by recalling the basic properties of the eigenmodes of the Dirac operator Ip 
and the hermitian operator 7^^ in the continuum and then discuss properties of the local 
chiral density which is the focus of this paper. Since {7^,^} = and Ip'^ = —Ip, the 
spectrum of Ip consists of pairs of eigenvectors ipx^x) and ip-xi^) = l^i^xi^) with imaginary 
eigenvalues ±iX. For A = the zero modes, {i'o^i}i<i<N, can also be chosen eigenstates of 
7^, 'y^ipo,i = iV'o.i since [7^, ^] ■j/'oi = 0. Here we have added the extra subscript i to the 
zero modes to allow for the possibility that more than one eigenvector with zero eigenvalue 
exists. The Atiyah-Singer index theorem requires that the number of zero modes, satisfy, 
N > l^l where u is the integer winding number of the gauge field (z^ = ... — 2, —1, 0, 1, 2, . . .). 

Clearly, zero modes of Ip are also zero modes of j^Ip. The remaining eigenvectors of 
7^(^ + m) can be easily related to those of ^ -|- m by diagonalizing the former on the 
two-dimensional subspace spanned by the basis {ipx,ilj-\}: 

^/TA \ f 0'\ ( m — i\ \ ( a\ 

-'(^ + '"'U = ™ + a U « 



The eigenvalues are easily seen to be 

\h = ± 

The normalized eigenvectors, iPh,Xh cire 




A^ + m^V (2) 



i'H.±\\H\ = — r li'm - i\) ^ i^x±i'm + i\) 2 ip_x\ . (3) 

V2(A2 + m2)4 L J 

Next we compare the chiral density as determined by these two sets of eigenvectors. For 
zero modes: 
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(4) 



the largest value possible, for both the hermitian and non-hermitian operators. For the 
other, paired eigenvectors and the non-hermitian operator, the identity ip-x{x) = 'y^'ipx{x) 
implies that product 

^{{xh^-^xix) = V^1;,(x)7^V^_a(x) = x{x) (5) 

is the same for each eigenvector in the pair. Similarly, the norm of the wave function 
at each point in space-time is the same for each pair of eigenvectors 

ip{{x)ipxix) = ipl^{x)i)^xix) = uj{x) (6) 

This is not true for the paired eigenvectors of the hermitian operator, where using Eq. |^ 
we find: 

'^H,XHi^)^H,\Hix) = w{x) ± ex{x) = ujh{x), (7) 

where 

^ (8) 

and A = ±|A|. In the chiral limit e vanishes, and the local chiral densities of the anti- 
hermitian and hermitian Dirac operators coincide. 

The above analysis should describe the properties of the low-lying eigenstates of the non- 
hermitian Wilson operator in the continuum limit. Thus, we conclude that the difference 
between the evaluation of the chiral density for the eigenvectors of the non-hermitian Wilson 
Dirac operator and those of the hermitian domain wall Dirac operator should agree in the 
continuum limit for those eigenvectors for which the parameter e is small. Thus, the results 
presented in this paper and those of Horvath, et al. address the same continuum question and 
can be compared. Since in our calculation the explicit fermion mass is very small (mj = 
for the Wilson action and nif = 5 ■ 10~^ for the Iwasaki action) we expect that for most 
modes the parameter e will also be very small making the normal and chiral densities the 
same for both operators. The only exception to this conclusion arises for the very smallest 
eigenvalues, A ~ 10~^ where e may be of order one due to the residual mass coming from 
mixing between the walls, rrij-es ~ 10~^ for the Wilson action and rrij-cs ~ 10~^ for the 
Iwasaki action P,pi][1.) 



III. DIAGONALIZATION METHOD 



Using the conjugate- gradient method proposed by Kalkreuter and Simma |T1], we cal- 
culate the 19 lowest eigenvalues and eigenvectors of Dh- This corresponds to the range 
^ A ^ 200 (Me V). In this section we describe our implementation of their methodQ. 



Code provided by Robert Edwards forms the basis for the computer program used here. 
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Following Kalkreuter and Simma, we use the conjugate gradient method to minimize the 
Ritz functional, 

on a sequence of subspaces orthogonal to the eigenvectors corresponding to the previously 
identified minima. When we have obtained our intended iVmax = 19 eigenvectors, we then 
evaluate the A^max x A^max matrix This small matrix is then diagonalized by 

a Jacobi method and the resulting transformation used to improve the initial set of vectors 
{^'fc}o<fc<Afmax- This combination of conjugate gradient minimizations and small matrix 
diagonalizations is repeated until a sufficiently accurate result has been obtained. 

This method requires the choice of two stopping criteria. For the sequence of A^max 
conjugate gradient minimizations of f^i'^k) we iterate each conjugate gradient procedure 
until the norm squared of the gradient of ni'^k) with respect to has decreased from 
its initial value by a factor of 10. This choice of 7 = 0.1, in the notation of Kalkreuter 
and Simma, was recommended by those authors and, after some testing, appeared to be a 
good choice for us. As a criterion to stop the outer iteration over A^max conjugate gradient 
minimizations followed by small matrix diagonalization, we required that the relative change 
in each of the eigenvalues be less than 10~^ after the final step in the method. 

Since we are actually interested in the small eigenvalues and corresponding eigenvectors 
of Dh, we performed a follow-on step after the above procedure has yielded a best set 
of eigenvectors and eigenvalues of Djj. In this step, we used the Jacobi method again to 
diagonalize the 19 x 19 matrix, and the resulting transformation to rotate the 

original basis of 19 eigenvectors into what should be eigenvectors of Dh- 

However, if the 19*^ and 20^^ eigenvalues of Dfj are nearly degenerate but the corre- 
sponding eigenvalues of Dh have opposite signs, then this 19- dimensional subspace will in 
general not be invariant under the application of Dh- The 19*^ eigenvector of Dj^ can be an 
arbitrary mixture of the 19*^ and 20*^ eigenvectors of Dh- If this happens, we expect to see 
one "spurious" and 18 valid eigenvectors of Dh after the Jacobi step. The spurious eigen- 
vector will be orthogonal to the other 18 and the magnitude of the corresponding eigenvalue 
will be distributed arbitrarily between a correct 19*^^ value and zero. 

One way to identify such a spurious eigenvalue is to look for an eigenvalue of Dh whose 
square does not agree with any of the 19 eigenvalues of Djj. In practice we order the eigenval- 
ues found for both Dh and Djj and consider 19 possible matches obtained by discarding one 
of the eigenvalues of Dh and comparing the squares of the remaining 18 with the lowest 18 
eigenvectors of D^. The case with the largest sum of squares of relative differences identifies 
the spurious eigenvalue. This algorithm produces 18 "good" eigenvalues and discards the 
largest eigenvalue in the case that no spurious eigenvector is present. We use this method 
to determine the 18 eigenvectors analyzed in this paper. 

A second method that in addition involves the eigenvector information evaluates the 
anti-commutator identity satisfied by Dh, 

{r5,D//} = 2myQ(-)+2g('^P), (10) 

where Q^"^^ and Q*-™^^ are pseudoscalar densities defined on the pairs of planes at s = 
{0,Ls — 1) and s = {Ls/2 — l,Ls/2), respectively, with precise definitions given in Ref. M. 
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We compute the diagonal elements of this matrix equation in our basis of eigenvectors of 

^n). (11) 

If for one or more eigenvectors the fractional difference of the left- and right-hand sides of 
Eq. |TT| is larger than 0.12 we discard the eigenvector producing the largest difference. If all 
eigenvectors give differences lying below this value, we discard the largest eigenvalue. The 
fraction 0.12 is a somewhat arbitrary dividing point that separates the few large discrepancies 
from the majority with a much smaller difference. 

Comparison of these two criteria shows disagreement on 4 of the 32 Wilson and 6 of the 55 
Iwasaki configurations. However, in each of those cases one method discards the nineteenth 
eigenvalue while the other identifies either the seventeenth or eighteenth as spurious. Since 
the eigenvectors of these largest eigenvalues are likely the least well determined, we view this 
as satisfactory and expect that the difficulties from this imprecise identification of "spurious" 
eigenvectors enter our sample of eigenvectors on the level of a few percent or less. 



)(mp) 



IV. ANALYSIS AND RESULTS 

In this section we examine some of the properties of the 18 lowest eigenvalues, Kn^i and 
the corresponding eigenvectors |A//i) determined by the procedure described above, 

DH\KH,i) = AH,i\^H,i) (0 < 2 < 18). (12) 

The connection between these eigenvectors and eigenvalues and those of the usual 4- 
dimensional Dirac operator was discussed in detail in 0| . The 5-dimensional wave functions, 
'^H,i{x, s) are expected to represent the wave functions of 4-dimensional Dirac eigenfunctions 
with an added exponential dependence on the fifth coordinate s, causing \'^H,i{x, s)p to fall 
rapidly as s moves away from the 4-dimensional planes s = and s = Lg — 1, as shown, for 
example in Figs. 29 and 30 of Ref Likewise the eigenvalues Ah,! should correspond to 
4-dimensional Dirac eigenvalues. 

For our conventions the components bound to the s = wall are predominately left- 
handed while those concentrated on the s = Lg — I wall are right-handed. It is important 
to recall that for a lattice with finite Lg, a left-handed mode can propagate with a very 
small probability to the right-handed wall, and visa-versa, giving a small breaking of chiral 
symmetry which, for low energy phenomena, can be described by a residual quark mass. 

Our analysis is based primarily on 55 quenched gauge configurations generated with the 
Iwasaki gauge action |]12[ at (3 = 2.6 and 4-dimensional lattice volume 16^. We also examine 



32 configurations generated with the standard Wilson action at /3 = 6.0 and the same 
volume. These two gauge actions have been choosen to correspond to approximately the 
same lattice spacing: ~ 2 GeV. The gauge configurations that we analyze are separated 
by 2000 sweeps of our updating algorithm, a simple two-subgroup heat-bath update of each 
link. The number of sites in the fifth dimension for the domain wall Dirac operator is 
Lg = 16, and the domain wall height is M5 = 1.8. For input quark mass we take rrif = 
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for the Wilson action and ruf = 5 ■ 10~^ for the Iwasaki case. This small, non-zero value of 
is used in the Iwasaki case to avoid the very slow convergence of the Kalkreuter Simma 
method that we found when using ruf = for that action. Finally, the residual quark mass 
for these couplings is quite small, mj-es ~ 0.0001 for Iwasaki P,|lO|, and rrires ~ 0.001 for 
Wilson which correspond to 0.4 and 4 MeV, respectively. 

To begin, let us examine the integrated or global chiral structure of Dh- We plot the 
magnitude of (Ajy^j|r^|Aj^j) in Figs. |I| and ^ for a typical "simple" and "complex" con- 
figuration. The matrix represents the physical 7^ matrix in the domain wall fermion 
formalism 

(AhA^'IAh,) = E ^L.(^' ^)sign(^ - (L. - l)/2)vl/^^,^(x, s) , (13) 

x£V s=0 

where V is the 4-dimensional space-time volume. 

The pattern seen in Fig. H for the simple configuration is precisely the chiral structure 
expected in the continuum. The single, diagonal element corresponding to \Ah,o) represents 
a zero mode which is an eigenstate of 7^ with eigenvalue +1. (We will refer to such modes 
as "near zero modes" since for our choice of parameters, their eigenvalues are not precisely 
zero.) In our entire sample of Iwasaki and Wilson configurations, all such near zero modes 
have very small eigenvalues, A ^ 10^'^ — 10~^, and all are either right-handed or left-handed 
within a given gauge configuration. Note this behavior is a natural consequence of the 
Atiyah-Singer theorem which requires an excess of right-handed to left-handed zero modes 
equal to the winding number of the background gauge configuration. This determines a 
minimum number of zero modes, all with the same chirality. The presence of additional 
zero modes would imply added constraints on the gauge background, corresponding to a 
set of zero measure if our near zero modes had a precisely zero eigenvalue. The remaining 
eigenvectors are grouped into pairs connected by 7^ precisely as expected for the continuum 
Dirac operator in the limit of vanishing mass. 

These simple continuum expectations are not satisfied by the complex configurations, 
such as that shown in Fig. |^. Of course, such configurations must be present for finite 
and finite lattice spacing. As the gauge configurations change continuously from one winding 
number to another, a plot of the sort shown in these figures must also change continuously 
and hence cannot always have the simple structure of Fig. While for the Wilson gauge 
action somewhat more than half of the configurations show the complex pattern of Fig. 0, for 
the better-behaved Iwasaki case, this fraction has dropped to 10%. It should be emphasized 
that all the low-lying eigenvectors studied, both complex and simple ones, fall off rapidly 
away from the walls with the minimum magnitude of the wave function between the walls 
falling at least a factor thirty below its value on the two physical boundaries. 

Such gauge configurations in which the winding number is changing can be associated 



with zero modes of the 4-dimensional Wilson Dirac operator with mass equal to — M5 ||13| , p!4 
Numerical simulations |]T5[ suggest that the density of such 4-dimensional Dirac operator 
zero modes decreases exponentially with the exponential of the coupling, p ~ e~'^^ so that 
such effects may vanish rapidly as the continuum limit is approached. 

Close to the continuum limit, the typical gauge configuration will be sufficiently contin- 



uous that its winding number can be identified [16]. As the winding number changes one 
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expects that localized, rapidly changing gauge fields will appear and small dislocations, on 
the scale of a very few lattice spacings will appear or disappear. It is natural to speculate 
that such configurations produce the complex matrix elements of Fig. ^ and the non-zero 
density of 4-dimensional Dirac zero modes described above. The comparison of Iwasaki and 
Wilson results suggests that while such configurations are quite common for the Wilson 
gauge action when ^ 2 GeV, they are dramatically suppressed under similar circum- 
stances by the form of the action proposed by Iwasaki. This general topic is the subject of 
much current research |[T7|-pO|. 

If we are to systematically evaluate the domain wall fermion, QCD path integral, we must 
include all configurations in our analysis. Although we are explicitly examining the small 
eigenmodes of the Dirac operator, the eigenvalues and eigenfunctions represent a complex 
mixture of both long-distance and short-distance physics. While we are able to explicitly 
focus on small Dirac eigenvalues, the background gauge fields contain the full spectrum of 
short- and long-distance fiuctuations. This is immediately demonstrated by the potential 
importance of small, very short-distance dislocations in the gauge configuration on quantities 
we are examining. However, we should also expect the more conventional short distance 
effects of wave function and mass renormalization, including short-distance contributions 
to the residual mass, to influence the low lying eigenvalues and eigenvectors analyzed here. 
A brief discussion of these effects on the eigenvalue spectrum can be found in Section VI 
C of Ref. [§]. We speculate that these complex configurations, which do not behave as is 
expected for smooth, continuum gauge fields, represent such short-distance effects. However, 
we should look carefully to see if configurations of the complex type introduce further chiral 
symmetry breaking at low energy beyond the simple residual mass described above. To date 
we have not recognized such effects. 

The global chirality of our eigenvectors is summarized in Fig. ^ The distribution shows 
a large narrow peak around zero corresponding to the non-zero modes and two smaller, but 
also narrow, peaks at ±1 corresponding to near zero modes. From these figures it is clear 
that we can easily distinguish near zero modes from non-zero modes, except for the handful 
of outliers with chirality neither close to zero nor ±1. 

Next we examine the distribution of local chirality, Xh{x)/Qh{x), which is shown in 
Figs. I and ^ for near zero modes and non-zero modes, respectively. Here the quantities 
Xh{x) and flnix) are the generalizations of the continuum quantities Xh{x) and ujh{x) 
defined in Eq. ^, to the case of domain wall fermions: 

0<s<Ls 

Xh{x)= sign(s-(L,-l)/2)vI/A^,(x,s)tv[/A^^(x,s) (14) 

0<s<Ls 

The six histograms superimposed in Figs. ^ and ^ correspond to histograms of the quantity 
^H{x)/flH{x) evaluated only at those sites where the normal density flnix) lies above a 
specified cut, for each of the eighteen eigenvectors that we have determined. The six cuts 
displayed correspond to the conditions: flnix) ■ 10^ > 3,4,5,6,7 and 8 which include 8.6, 
4.3, 2.4, 1.5 1.0 and 0.7% of the sites in the space-time lattice, respectively. These cuts 
correspond to sampling the local chirality, on average, from sites that account for 28, 19, 13, 
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10, 8, and 6%, respectively, of the total probability density of each eigenvector. In Table | 
we give these numbers with more precision for both the Iwasaki and Wilson cases. 

The near zero mode distributions are sharply peaked at ±1 for all cuts, as expected. 
The non-zero mode distributions are also clearly double-peaked, with peaks centered ap- 
proximately around ±.95 —>■ ±.8, depending on the size of the cut. The distributions fill in 
between the peaks as more sites are sampled. Thus, our data show that definite chirality is 
strongly correlated with local maxima of both near zero and non-zero eigenmodes. The latter 
is certainly consistent with an instanton-dominated vacuum picture and is in disagreement 
with the recent work of Horvath, et ai. 

We also note that the strong correlation between chiral and normal density for our non- 
zero mode eigenvectors is also evident if we use the definition of local chirality in instead 
of the ratio Xh{x)/Qh{x). In Fig. ^ we show a similar set of histograms for the Wilson 
gauge action. Clearly the same strong correlation between chiral and normal density is seen 
for this case as well. 

We now discuss two important consistency checks on this potentially interesting result. 
First we determine the distribution of "global" chirality of each eigenvector when computed 
by including only those lattice sites selected by our cuts on flnix). We expect this distribu- 
tion to be similar to that found when summing over all sites in the lattice, with each non-zero 
mode showing approximately zero total chirality. The results are shown in Fig. 0. While 
considerably broader than the distribution seen for the global chirality, it is still strongly 
peaked around zero with the sharpness of the peak decreasing as the cut on Qh{x) is made 
more stringent. This broadening is easily understood as the statistical effect of examining 
the sum over a smaller number of space-time points. Over-all, Fig. |^ is quite consistent 
with the expectation that chirality of the non-zero mode eigenvectors is quite evenly split 
between left- and right-handed lumps. 

As a second consistency check, we investigate the extent to which the lowest lying 18 
eigenmodes which we examine actually play an important role in the low energy physics 
described by the gauge configurations being studied. An easy way to address this question is 
to compute the contribution of the modes which we have isolated to the chiral condensate and 
to compare that contribution to the total chiral condensate determined independently [Q. 
In Table |I| and Fig. ^ we present such a comparison. As can be seen, these lowest 18 
eigenvectors provide a large fraction of the actual value of {qq) for the light input mass, 
TUf = 0. For larger values of ruf, the contribution of these modes is a rapidly decreasing 
fraction of {qq) . However, for non-zero quark mass, the quadratically divergent contribution 
of high mass states oc mf/a?' should make an increasingly important contribution. We can 
easily include these effects, by using a fit to (gg) of the form 

- = ^ hao + aim/, (15) 

where a_i, Oq, Oi and dmi^qq) are parameters determined and tabulated in Ref. @]. The 
large coefficient ai describes these divergent contributions, allowing them to be included in 
our estimate of (gg) by simply adding the term ai ■ (m^ + m^es) to —{qq). These results 
are included in Table and Fig. |^. One sees that when this expected contribution from 
far off-shell states is included, we have good agreement with the directly computed values 
(gg) for m/ in the range between 0.0 and 0.01. This suggests that the 18 lowest modes that 



9 



we have examined provide the bulk of the physical vacuum expectation value (gg) and are 
therefore relevant to an attempt to characterize the physics of the QCD vacuum. 

Finally, we note that the above results for the chiral density distribution of the near zero 
modes shown in Fig. | are not very illuminating since the local chirality is either left- or 
right-handed for all points in the lattice for these eigenvectors. A more informative picture 
is given in Fig. ^ which shows the magnitude of the right-handed components of a near 
zero mode in the x — y plane summed over z and t. The mode is clearly localized (and 
remains so in the other planes that are not shown). It is also dominantly right-handed since 
the left-handed component shown in Fig. |TU| is more than 50 times smaller. These figures 
corresponds to the unique, near zero mode labeled as eigenvector in Fig. |I[ 

For completeness, we include similar figures showing the spatial distribution of the right- 
(Fig |1^) and left-handed (Fig. |T2|) parts of the first non-zero mode, numbered 1 in Fig. |1|. 
These provide a particular illustration of the behavior shown in average terms in Fig. |^. 
The right-handed components of eigenvector 1 are centered around the point (x, t) = 
(13, 0, 7, 15) while the left-handed components of this eigenvector are largest at the distinct 
point (11,9,3,6). This specific case shows a clear correlation between localization and 
chirality as would be expected in an instanton picture. 



V. CONCLUSION 

We have examined the local correlation between the chirality and magnitude of the low- 
lying eigenmodes of the Dirac operator in quenched lattice QCD. A strong correlation is 
observed in our calculation suggesting a picture of the QCD vacuum containing background 
gauge fluctuations producing space-time regions in which the Dirac eigenstates are both 
localized and chiral. For those 0.7% of the lattice sites at which 6% of the eigenmode mag- 
nitude is concentrated, these eigenmodes are nearly eigenstates of 7^ with XH{x)/flH{x) 
within 90% of 1. While this is certainly consistent with a vacuum described as a superpo- 
sition of instantons, it is an important feature that may be accommodated by other models 
as well. 

An essential ingredient of our calculation is the use of the domain wall fermion formu- 
lation. The domain wall Dirac operator readily identifies topological zero modes on gauge 
configurations as they appear in the Feynman path integral without resort to smoothing 
or cooling. The presence of these near zero modes in the domain wall Dirac spectrum is 
easily seen in the small mass behavior of the chiral condensate and their implications for 



quenched calculations of the both the chiral condensate and hadron masses was discussed in 
detail in Ref. |Q. Related behavior has been found using the overlap Dirac operator P^ , p5 |, 
an alternative lattice fermion formulation, also with improved chiral symmetry. 

The recent calculation [jl| using Wilson fermions which motivated this study found no 
chirality-localization correlation. However, the calculation of the present paper is improved 
over this earlier work in two significant respects. The most important is that we use domain 
wall fermions which preserve chiral symmetry to a high degree at non-zero lattice spac- 
ing while Wilson fermions do not. In addition the lattice spacing used here is roughly 
twice as small as that of Ref. |I| . 

Note added: After this paper was essentially complete, the recent work of DeGrand and 
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Hasenfratz became available p4[. Our analysis and conclusions are quite similar to those 
of this paper. While both that paper and the present one, use improved fermion Dirac 
operators, overlap and domain wall respectively, we have used the gauge configurations 
directly, without smearing or fattening. This suggests that this strong correlation between 
spacial localization and chirality is seen even for gauge configurations whose fluctuations span 
the full range of distance scales produced by both the Iwasaki and Wilson gauge actions. 
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TABLES 

TABLE L Fractions of lattice sites and eigenvector normalization included for the six cuts, 
> rimiii) used in this paper. 





Iwasaki Action 


Wilson Action 


Q ■ • 1 


Fraction of 


Fraction of 


Fraction of 


Fraction of 




lattice sites 


normalization 


lattice sites 


normalization 


3 


0.086 


0.285 


0.080 


0.319 


4 


0.043 


0.186 


0.042 


0.235 


5 


0.024 


0.132 


0.026 


0.186 


6 


0.015 


0.100 


0.017 


0.156 


7 


0.010 


0.079 


0.012 


0.136 


8 


0.007 


0.064 


0.009 


0.121 



TABLE II. Contributions to the chiral condensate from the 18 lowest eigenvectors examined 
in this paper, {qq)is^ for five different quark mass values. In addition, we include the ~ 1/a^ 
contribution from the large eigenmodes by including the term ai(m/ + mres) using the values 
ai = 3.97(8)10"^ and rrires = 1.24(5)10"^ taken from Ref. Q. Finally these results are compared 
with the complete value for (qq) found in that paper. This comparison is done in the case of the 
Wilson action at /3 = 6.0 where the detailed results for (qq) are available from that earlier work. 



rrif -{qq) 18 -(gg)i8 + ai{mf + mres) -{qq) 

0.0 0.0036(6) 0.0037(6) 0.00219(20) 

0.0025 0.00087(9) 0.00119(9) 

0.005 0.00059(5) 0.00113(5) 0.00100(2) 

0.0075 0.00048(3) 0.00124(3) 

0.01 0.00043(2) 0.00140(2) 0.00134(1) 
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FIGURES 




FIG. 1. The magnitude of the matrix elements {AH,i\^^\^H,j) evaluated on a "simple" config- 
uration. The pattern seen is close to that expected for matrix elements of 7^ in the continuum 
theory: The single near zero mode is an an approximate eigenstate of F^ while the non-zero modes 
come in pairs related by F^^Ah = "^A^- 
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FIG. 2. The magnitude of the matrix elements {AH,i\T^\^H,j) evaluated on a "complex" con- 
figuration. For the Wilson gauge action more than one-half of the 32 configurations look like this 
while for the Iwasaki action only one in ten has this complex structure with the remaining 90% 
appearing similar to Fig. ||. 
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FIG. 3. The distribution of global chirality of the eigenvectors of Dh evaluated at zero input 
quark mass, mj = 0. 
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FIG. 4. The distribution of local chirality of near zero mode eigenvectors of Dh 

on the Iwasaki ensemble of gauge fields for sites where ^^^(x)^(x) is greater than an arbitrary 
imposed cut. The different cuts correspond to keeping between 1 and 8% of the total sites in the 
space-time lattice. 
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FIG. 5. The same quantities as in Fig. ^, but for non-zero mode eigenvectors and again for the 
Iwasaki action. The double-peak structure is a feature expected in instanton-dominated models of 
the QCD vacuum. 
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FIG. 6. The same quantities as in Fig. ^, except now the Wilson gauge configurations are 
examined. The double-peak structure is very similar to that found in the Iwasaki case. 
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FIG. 7. The distribution of chiral density for non-zero mode eigenvectors of Dh summed over 
those lattice sites obeying the cuts imposed on ^Ih{x) as described in the text. As the cut is 
increased, the width of the distribution becomes more narrow. The distributions indicate that on 
average, chirahty for a single eigenvector is evenly split into left- and right-handed lumps. 
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FIG. 8. (qq) for quenched simulations done on 16^ X 32 lattices at /3 = 6.0 for = 16. In 
addition to the direct results for (qq) presented in Ref . Q , and the fit of the form given in Eq. 15 , 
we also plot the values obtained from the 18 eigenvectors studied here (filled triangles) and the 
sum of those values with the contribution from far off-shell states, ai(mj + mros) (filled diamonds). 
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FIG. 9. The distribution of the magnitude of the right-handed components of a near zero mode 
eigenvector of Dh in the x — y plane, summed over the z and t coordinates. A single, well-localized 
peak is visible near a; = 14 and y = 2 — 3. 
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FIG. 10. The distribution of the magnitude of the left-handed components of a near zero mode 
eigenvector of Dh in the x — y plane, summed over the z and t coordinates. This is the same 
eigenvector whose right-handed components are show in Fig. While a peak can also be seen in 
this figure, its magnitude is 50 times smaller than that shown for the right-handed components. 
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FIG. 11. The distribution of the magnitude of the right-handed components of a non-zero 
eigenvector of Dh in the x — y plane, summed over the z and t coordinates. A single somewhat 
broad peak is seen with x = 13 and y = or 15. [Note: it should be remembered that our periodic 
boundary conditions imply that the points y = and y = 15 are adjacent.] 
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FIG. 12. The distribution of the left-handed components of the non-zero eigenvector shown in 
Fig. 11. Note this shows a quite different peak from that in the previous figure with (x, y) = (11, 11) 
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